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ABSTRACT 

An analysis of high-resolution near-infrared spectra of a sample of 45 asymptotic giant branch 
(AGB) stars towards the Galactic bulge is presented. The sample consists of two subsamples, 
a larger one in the inner and intermediate bulge, and a smaller one in the outer bulge. The 
data are analysed with the help of hydrostatic model atmospheres and spectral synthesis. 
We derive the radial velocity of all stars, and the atmospheric chemical mix ([Fe/H], C/O, 
Al, Si, Ti, and Y) where possible. Our ability to model the spectra is mainly limited 
by the (in)completeness of atomic and molecular line lists, at least for temperatures down 
to Teff ~ 3100K. We find that the subsample in the inner and intermediate bulge is quite 
homogeneous, with a slightly sub-solar mean metallicity and only few stars with super-solar 
metallicity, in agreement with previous studies of non-variable M-type giants in the bulge. 
All sample stars are oxygen-rich, C/O<1.0. The C/O and carbon isotopic ratios suggest that 
third dredge-up (3DUP) is absent among the sample stars, except for two stars in the outer 
bulge that are known to contain technetium. These stars are also more metal-poor than the 
stars in the intermediate or inner bulge. Current stellar masses are determined from linear 
pulsation models. The masses, metallicities and 3DUP behaviour are compared to AGB evo¬ 
lutionary models. We conclude that these models are partly in conflict with our observations. 
Furthermore, we conclude that the stars in the inner and intermediate bulge belong to a more 
metal-rich population that follows bar-like kinematics, whereas the stars in the outer bulge 
belong to the metal-poor, spheroidal bulge population. 

Key words: stars; AGB and post-AGB - stars; late-type - galaxy; bulge. 


1 INTRODUCTION 


The Galactic bulge is a complex stmcture in the central regions of 
the Milky Way Galaxy. Its stellar populations have been studied 
in some detail in the recent years. Stellar types that have been in¬ 
vestigated by high-resolution spectrosc opy in the past range from 


microlensed dwarf and subgiant stars (lBensb^_e^alJj201^^_and 
references therein), to K-type giants (e .g. Babusiam(_et_al.l I 20 1 d ; 
iHill et al.ll2^; lUttenthaler et alj|2012l; iNess et alj|2013h . to M- 
type giants dRich et alJl2012L and references therein). Also plan- 
etary nebulae (PNe) in the bulge have been studied spectroscopi- 


* Stefan.uttenthaler@univie.ac.at (SU); joris.blommaert@vub.ac.be 
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cally jMaciell[19991: IChiappini et al.|[2Q0^ : iGuzman-Ramirez et al.l 

12011 . 

A group of bulge stars that has so far not been studied in de¬ 
tail by high-resolution spectroscopy consist of the asymptotic giant 
branch (AGB) stars. AGB stars represent the last stage of stellar 
evolution where the energy output is dominated by nuclear burning. 
This stage is characterised by strong pulsations of the envelope that 
enhance the stellar wind, which is accelerated by radiation pressure 
acting on dust particles formed in the outflowing material. Eventu¬ 
ally, the stellar evolution is terminated by the high mass loss and 
the star may enter the PN phase. In the deep interior of the star, 
nucleosynthetic processes take place that produce carbon ('^C) and 
heavy elements. These may subsequently be brought to the sur¬ 
face by a deep mixing event called the third dredge-up (3DUP). If 
enough C is added to the atmosphere, the abundance ratio of C to 
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O (by number of atoms) will exceed 1.0 and the star will become 
a carbon star (C-star). The formation of C-stars depends on the ini¬ 
tial mass and metallicity; A star may turn into a C-star only if it is 
massive enough to undergo enough 3DUP events, where the lower 
mass limit is a function of metallicity, be ing lower for lower metal¬ 
licity jMarigo et al.l2013l:lKaraka^ 20]jb. The lack of C -stars in the 

20131) suggests that there could 


■ regions of M31 Isover et a~ 


be a metallicity threshold above which the formation of C-stars is 
inhibited altogether. 

AGB stars in the Galactic bulge are interesting for a number 
of reasons. First of all, it has been shown that intrinsic carbon stars, 
i.e. stars that owe their enhancement of C to internal nucleosynthe¬ 
sis and 3DUP and not to binary mass transfer, are widely absent 
in the Galactic bulg e ( Azzo£Mdi_et^ 19961: ^anco & Temdrud 


ll989l:lTvson & Riclj|l99ll : lNg||l997l:ISchultheislll998h . This would 

mean that either the mass of bulge stars is too low or the metallicity 
too high to undergo (sufficient) 3DUP to form C-stars. However, it 
is known that oxygen-rich AGB stars in the outer Galactic bulge 
show absorption lines of technetium (Tc), a c lear indicator of re¬ 
cent or ongoing 3DUP dUttenthaler et aklEoPTh . This suggests that 
3DUP does happen in bulge stars, but probably not enough to form 
C-stars. An exception may be the first carbon Mira in the direction 
of the Galactic bulge that was recently presented bv iMiszalski et al.l 
ll2013h . It is not clear if that star is indeed a member of the Galactic 
bulgfl 

Samples of non-variable M-type gi ants, the precursors o f 
variable AGB sta rs, h ave been studied by iRich & Origli^ l l2005h . 
iRich et af] ( l2007h . and iRich et al.l ( l2012h . One of their main results 
is that the mean metallicity of M-giants in the bulge is slightly 
sub-solar, and that stars with super-solar metall icity a r e not com- 
mon. On the other hand, it was concluded by iBue i 1 I2OI3I) that 
PNe in the Galactic bulge, the evolutionary successors of AGB 
stars, have a re latively high, i.e. on ave rage super-solar, metallic¬ 
ity. Yet again, lUttenthaler et alj | |2012|) suggests that AGB stars 
in the outer bulge predominantly descend from the metal-poor 
population ([Fe/H] » -0.6), not from the metal-rich population 
([Fe/H] K -1-0.3) of the bulge. This latter result was derived from 
the fact that the radial velocity dispersion of AGB stars and also of 
PNe in the outer bulge is similar to that of the metal-poor popula¬ 
tion, which differs significantly from that of the metal-rich popula¬ 
tion. _ _ 

Finally, _ IPerea-Calderon et ^ ( l2009h and 

iGuzman-Ramirez et al. ( I2OIII) showed that many Galactic 
bulge PNe exhibit evidence of mixed chemistry with emission 
from both silicate dust and polycyclic aromatic hydrocarbons 
(PAHs), the latter being a signature of C-rich chemistry. Such 
signatures are difficult to explain in an environment with ( assum ed) 
low C/0 ratio, which required iGuzman-Ramirez et al.l ( 1201 ih to 
invoke hydrocarbon chemistry in an UV-irradiated, dense torus 
around the central stars of the PNe. In their chemical model, the 
formation of PAHs would be considerably eased if the C/O ratio 
was increased from its initial value by at least a few 3DUP events. 

These considerations make Galactic bulge AGB stars interest¬ 
ing targets for further investigation. A determination of the metal¬ 
licity and mass of AGB stars would help in understanding the lack 
of intrinsic C-stars in the bulge. In particular, a metallicity determi- 


* It was reported bv ICole & Weinberg l2002h that the Galactic bar is traced 
by luminous, very red (in J — K) stars. These were suspected by Cole & 
Weinberg to be infrared carbon stars. However, these objects lack spectro¬ 
scopic confirmation of their C-star nature. 


nation would shed light on the discrepancies in the metallicities of 
non-variable M-type giants and PNe in the bulge. A measurement 
of the C/O ratio in bulge AGB stars would provide observational 
constraints to explanations for the mixed-chemistry phenomenon 
in Galactic bulge PNe. 

In this paper we present near-IR spectroscopic observations 
of such a sample of stars and the chemical composition derived 
from them. The near-IR range has a number of advantages over the 
optical spectral range for studies of bulge red giants, most of all 
the lower interstellar extinction, less crowding of lines, and smaller 
effects of departu res from local thermodynamic equilibrium (see 
iRvde et alJiBol for a detailed discussion). Als o 3D effects are 
small er in the near-IR than in the optical range jKucinskas et al.l 
l20I3h . Thus, with these data, abundance determinations are done 
that w ould not be possible to obtain from optical spectra, as used 
e.g. in lUttenthaler et af] ( l2007h . 

The paper is organised as follows: The sample selection, ob¬ 
servations, and data reduction are presented in Sect.|^ the analysis 
of that data is described in Sect.[3j the results and their discussion 
are to be found in Sect.|4j and conclusions are drawn in Sect.|5] 


2 SAMPLE SELECTION, OBSERVATIONS, AND DATA 
REDUCTION 

In order to study the dust formation in the circumstellar environ¬ 
ment of AGB star s, a mid-infrared sp ectroscopic survey with the 
IRS spectrograph dHouck et al.ll2004) on board the Spitzer space 
telescope was performed on a sample of Galactic bulge AGB stars 
(Spitzer Proposal ID #3167, PI J. Blomma ert; iBlommaert et all 
I2OO7I : I^nhollebek3l2007l : IgoMz et alll20I4l) . The targets for tha t 
survey were selected from the ISOGAL survey dOmont et alj2003h 
in intermediate bulge fields (|fi| = 1° to 4°, corresponding distance 
from the mid-plane ~ 140-560 pc at Ro = 8.0 kpc). The (As-[15]) 
colour, where [15] is the 15^m hand of the IS O satellite, is an 
estimate of the dust m ass-loss rate of the stars dOiha et '^ 120031: 
IBlommaert et alj[2006h and was used to ensure that the full range 
of mass-loss rates, starting from the onset of (dust) mass loss up 
to the mass-loss rates of the so-called super-winds of OH/TR stars 
(10“‘*MQ/yr) was included in the sample. 

We selected from this Spitzer IRS sample those stars for high- 
resolution near-IR observations for which we expected to he able to 
model their near-IR spectra, i.e. OH/IR stars were excluded. Nev¬ 
ertheless, four large-amplitude va riables (LAVs) were ob served. 
Following the recommendations of lMcSavenev et akl d2007h . these 
four LAVs were selected such that they would he near minimum 
light at the time of the observing run. The initial goal of the ob¬ 
servations was to determine the C/O and '^C/'^C ratios of the 
stars as an “evolutionary clock” along the AGB as a result of re¬ 
peated 3DUP. The observations were carried out in visitor mode 
(observer: J. Blommaert) with the CRy ogenic high-resol ution In- 
fraRed Echelle Spectrograph (CRIRES iKaufl et ^l2004ll in three 
nights between 09 and 12 June 2008. The target stars themselves 
were used as guide star for the wave-front sensor of the adaptive 
optics system. An entrance slit width of 0.4"was chosen, which re¬ 
sults in a nominal resolving power of R = A/AA = 50 000. 

Along with the 37 bulge AGB stars from the Spitzer sample, 
spectra were also obtained of a few comparison stars that had been 
studied previously with high-resolution near-IR spectroscopy by 
other authors. These include BMB 78 stu died by ICunha & SmithI 
(l2006l) . BD-0I297I, which was studied bv iRich & Origlial 1 2005 1. 
as well as BMB 289, which was investigated by both these works. 
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Furthermore, spectra of five K/M-type field giants were obtained; 
the analysis and results for the field stars are reported in the Ap¬ 
pendix. 

Also included in the present paper is an analysis of CRIRES 
spect ra of eight A GB stars in the outer bulge, the so-called Plaut 
stars ([pi 3U97II). They are much further from the Galactic mid¬ 
plane than the Spitzer stars, at ~ -10° or ~ 1.4kpc. These 
stars w ere selected from the sample discussed in lUttenthaler et al] 
(l2007h and include all four stars that were found to contain Tc in 
their atmosphere, as well as four Tc-poor stars. These observations 
were done either in the science verification runs (August 2006) or 
in a regular service mode observing programme (Prog. ID 383.D- 
0685(A), September 2009). 

All stars were observed at least in one setting in the H-band 
(grating order 36) and in the K-band (grating order 24). The H-band 
spectra cover approximately the range Tyac - 1543 to 1561 nnfl 
with a gap of ~ 2 nm between the detectors. Due to contamination 
from neighbouring orders, only data from detector chips 2 and 3 are 
considered here. The H-band spectra are dominated by the '^CO 
3-0 band head, a number of OH and CN lines, as well as atomic 
lines (Fe, Ni, Si, Ti, S). The K-band spectra in grating order 24 
cover the range dvac - 2355 to 2407 nm and are dominated by CO 
Av = 2 lines. 

Sixteen of the science targets and one of the comparison stars 
were observed in two additional settings in the K-band to also have 
lines available of atomic species that are relevant for the dust for¬ 
mation. These spectra in the grating orders 27 and 26 cover the 
range Tvaj, 2087 to 2134nm and 2174 - 2222 nm, respectively. 
Besides being cluttered with abundant CN lines, these spectra con¬ 
tain, amongst others, lines of Na, Mg, Al, Si, Sc, Ti, and Fe. The Na 
lines were targeted to check a possible correlation between Na line 
equivalent width measured from low-resolution spectra and dust 
mineralogy in the circumstellar envelope. This correlation was not 
reproduced with the high-resolution spectra. 

A basic reduction of the data was done with the standard ESC 
CRIRES pipeline, version 1.7.0. Eurther processing, in particu¬ 
lar wavelength calibration and telluric correction, was done with 
custom-made IDL tools. The H-band spectra were wavelength- 
calibrated manually using the ThAr lamp exposures taken at day 
time. Due to the lack of suitable ThAr lines in the K-band settings, 
these spectra were wavelength-calibrated using the numerous tel¬ 
luric features imprinted in them. 

Before dividing the science target spectrum by the telluric 
standard star spectrum for telluric correction, both spectra were 
re-binned to a common wavelength vector and the airmass of the 
telluric spectrum was adjusted to that of the corresponding science 
spectrum according to the Lambert-Beer law: 


./std.adj — exp 


ln(/s,d)^ 

-^std 


with /sid the flux of the telluric standard star spectrum continuum- 
normalised to 1.0, and Xsa and Astd the airmass of the science target 
and of the telluric standard star, respectively. The science spectrum 
was then divided by the adjusted standard star spectrum. This care¬ 
ful procedure in the telluric correction yielded very satisfactory re¬ 
sults. Weak residuals from the division were only discernible in 
the cores of the strongest telluric lines. No telluric correction was 
applied to the H-band spectra because the wavelength range un¬ 
der consideration is essentially free from telluric absorption. The 


achieved signal-to-noise ratio (S/N) per pixel of the reduced one¬ 
dimensional spectra is relatively high (~ 100) for most stars, only 
for a few stars it is somewhat low (of the order of 30). The S/N is 
thus not the limiting factor in the analysis of most stars. 

Particular care was taken in the final step of data reduction, the 
continuum normalisation. The prime method used in this step was 
to identify narrow pieces of spectrum (a few 0.01 up to O.lOnm) 
that are basically unaffected by line absorption in as many of the 
sample stars as possible. The mean flux of the 30% brightest pix¬ 
els in these ranges was set to define the continuum (flux=1.0). This 
continuum placement was always checked in comparison with syn¬ 
thetic spectra. Eigure|5] provides an example of how well the con¬ 
tinuum in observed and synthetic spectra are matched. In some 
cases, in particular for the coolest and most strongly variable stars, 
even these narrow pieces of spectrum were affected by line absorp¬ 
tion. In these cases, the continuum was placed manually, mostly by 
dividing it by its maximum flux, except in cases that were obvi¬ 
ously affected by bad pixels, residuals of the telluric correction, or 
emission components. This method was necessary anyway in the 
K-band spectra of grating order 24 because this spectral region is 
heavily cluttered with CO Av = 2 lines. 


3 ANALYSIS 

3.1 Stellar parameter determination 

Before determining elemental abundances from the high-resolution 
spectra using spectral synthesis techniques, the main physical pa¬ 
rameters of our sample stars need to be determined. We derived 
those stellar parameters from low-resolution optical spectra, broad 
band photometry, and the pulsation periods of the stars. These data 
will be described in more detail in a separate paper (Vanhollebeke 
et al., in preparation). 


3.1.1 Pulsation period 

The pulsation periods of the sample stars are an important ingre¬ 
dient for the stellar parameter determination. K-band photometry 
of the Spitzer sub-sample span ning ~ 820 days wa s obtained by 
one of us (PW) and analysed bv IVanhollebeka ( l2007h . who also de¬ 
rived pulsation periods where possible. Eurthermore, the VizieB0 
database was searched for additional period determinations. The 
Optical Gravitational Lensing Experiment (OGLE)-II catalogue 
provided pulsation periods for a few more objects, in particular 
among the short-period, low-amplitude variables. In cases w here 
two periods were available, the one from IVanhollebekj ( l2007h was 
adopted. Eor the Plaut stars, the pe r iods a nd K- band magnitudes 
were taken from lUttenthaler et al.l l l2007h and lUttenthaler et al.l 
( l2008h . The mean K-band magnitudes and adopted periods are 
listed in columns 4 and 5, respectively, of Table [T] 


3.1.2 Temperature 

The most important stellar parameter for our spectroscopic anal¬ 
ysis is the effective temperature. To this end, low-resolution flux- 
calibrated optical spectra were obtained by one of us (PW) with 
the Double-Beam Spectrograph (DBS) mounted to the Australian 
National University 2.3 m telescope at Siding Spring Observatory 


^ All wavelengths in this paper are expressed as nm in vacuum. 


^ http://vizier.u-strasbg.£r/viz-bin/VizieR 
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Table 1. Main parameters of the sample stars. Meaning of the columns: object name; Galactic longitude /; Galactic latitude b\ mean 
K-band magnitude; pulsation period; adopted effective temperature; luminosity der ived by integrating the SED; luminosity con'ected 
for bulge depth scatter using the period - K-magnitude relations of llta et ^ j20Q4l see Sect. |3.1.3t ; current pulsation mass; adopted 
log g of the model atmosphere used in the analysis (see Sect. [HI; LAV: large-amplitude variable. Unknown values are replaced by 
The luminosities of the comparison stai's BMB 78, BMB 289, and BD—012971 are calculated from their temperatures and surface 
gravities, assuming 1 Mq. 


Object name^ 

1 

[deg] 

b 

[deg] 

{K) 

(mag) 

Period 

(days) 

TeB 

(K) 

L^/Lq 

uncorr. 

LJLq 

corr. 

M/M© 

cgs 

Remark 

BMB 78“ 

1.15 

-3.78 

7.53 


3600 

654 



0.80 


BMB 289“ 

1.22 

-4.02 

6.17 


3375 

1268 



0.40 


BMB 289’’ 



6.17 


3200 

814 



0.50 


BD-012971'’ 

-11.5 

50.8 

4.29 


3600 

1304 



0.50 


J174117.5-282957 

-0.13 

1.04 

7.10 


2800 

5609 



-0.28 


J174123.6-282723 

-0.08 

1.04 

8.32 


3200 

1638 



0.13 


J174127.3-282851 

-0.09 

1.02 

7.26 

438.95 

2600 

4732 


1.85 

-0.26 

LAV 

J174127.9-282816 

-0.08 

1.02 

7.14 


3000 

4921 



-0.20 


J174128.5-282733 

-0.07 

1.02 

7.48 

309.46 

3000 

3652 

4804 

1.37 

-0.20 


J174139.5-282428 

-0.01 

1.02 

7.04 


2800 

5314 



-0.28 


J174140.0-282521 

-0.02 

1.01 

7.26 


2900 

4813 



-0.20 


J174155.3-281638 

0.14 

1.04 

6.82 


2700 

6324 



-0.28 


J174157.6-282237 

0.06 

0.98 

7.58 


3300 

3342 



-0.13 


J174203.7-281729 

0.14 

1.00 

7.55 

392.12 

2800 

4541 


1.38 

-0.23 

LAV 

J174206.9-281832 

0.13 

0.98 

7.26 

481.61 

2800 

5152 


1.29 

-0.28 

LAV 

J174917.0-293502 

-0.14 

-1.02 

8.55 


3200 

1541 



0.16 


J174924.1-293522 

-0.13 

-1.04 

8.54 


3700 

1680 



0.37 


J174943.7-292154 

0.10 

-0.99 

7.94 


3000 

2413 



-0.15 


J174951.7-292108 

0.12 

-1.01 

8.26 


3500 

1795 



0.25 


J175432.0-295326 

0.18 

-2.16 

7.29 

401.74 

2900 

3844 


0.94 

-0.20 


J175456.8-294157 

0.39 

-2.14 

6.90 

322.11 

2900 

5699 

5410 

1.82 

-0.23 


J175459.0-294701 

0.32 

-2.19 

6.96 

472.15 

2600 

6095 


2.38 

-0.26 

LAV 

J175515.4-294122 

0.43 

-2.19 

8.26 

76.90 

3300 

1652 

2059 

0.69 

0.00 


J175517.0-294131 

0.43 

-2.20 

7.99 


3500 

2483 



0.11 


J180238.8-295954 

0.96 

-3.73 

7.86 

119.00 

3200 

2005 

3843 

0.92 

-0.21 


J180248.9-295430 

1.05 

-3.72 

8.13 

75.69 

3300 

1759 

2116 

0.68 

0.00 


J 180249.5-295853 

0.99 

-3.76 

8.65 

65.18 

3400 

1214 

1848 

0.73 

0.12 


J180259.6-300254 

0.95 

-3.82 

7.46 

228.87 

3000 

2971 

3374 

1.22 

-0.14 


J180301.6-300001 

1.00 

-3.81 

8.37 

33.51 

3400 

1330 

606 

0.45 

0.38 


J180304.8-295258 

1.10 

-3.76 

8.53 

35.52 

3400 

1147 

683 

0.48 

0.38 


J180305.3-295515 

1.07 

-3.78 

7.15 

207.02 

3100 

4023 

2944 

1.02 

-0.10 


J180305.4-295527 

1.07 

-3.78 

8.00 

113.94 

3200 

2015 

4184 

1.15 

-0.14 


J180308.2-295747 

1.04 

-3.81 

7.92 

113.04 

2800 

1919 

3613 

2.13 

-0.04 


J180308.7-295220 

1.12 

-3.77 

7.94 

98.00 

3100 

1782 

2842 

1.04 

-0.07 


J180311.5-295747 

1.05 

-3.82 

7.39 

292.04 

2900 

3036 

4557 

1.64 

-0.20 


J180313.9-295621 

1.07 

-3.82 

8.04 

75.00 

3400 

1801 

2073 

0.61 

0.05 


J180316.1-295538 

1.09 

-3.82 

8.31 

62.69 

3300 

1328 

1519 

0.65 

0.13 


J180323.9-295410 

1.12 

-3.83 

8.15 

84.84 

3200 

1547 

2409 

0.82 

-0.05 


J180328.4-295545 

1.11 

-3.86 

7.40 

119.09 

3200 

3005 

3850 

1.15 

-0.11 


J180333.3-295911 

1.06 

-3.90 

8.44 

203.95 

3300 

1247 

2970 

0.78 

-0.11 


J180334.1-295958 

1.05 

-3.91 

7.06 

119.09 

3100 

3986 

3696 

1.06 

-0.18 


Plaut 3-45 

-0.94 

-7.22 

6.65 

271.02 

3000 

5308 

3938 

1.25 

-0.19 

LAV 

Plaut 3-100 

-3.85 

-8.92 

6.37 

298.70 

3000 

8288 

5432 

1.68 

-0.21 

LAV 

Plaut3-315 

0.04 

-7.65 

6.54 

326.80 

3100 

6351 

5514 

1.30 

-0.27 

LAV 

Plaut 3-626 

-2.63 

-10.10 

7.19 

298.48 

3400 

3835 

5344 

0.89 

-0.26 

LAV Tc-rich 

Plaut 3-794 

0.92 

-9.03 

6.04 

303.54 

3000 

11278 

5584 

1.71 

-0.21 

LAV 

Plaut 3-942 

-2.41 

-11.47 

6.55 

338.00 

3500 

6753 

6235 

0.83 

-0.31 

LAV Tc-rich 

Plaut 3-1147 

-1.84 

-12.31 

5.77 

395.63 

3000 

13506 

7577 

1.84 

-0.31 

LAV Tc-rich 

Plaut3-1347 

0.35 

-12.60 

5.99 

426.60 

3300 

10730 

8143 

1.16 

-0.38 

LAV Tc-rich 


^ Stellar pai'ameters from ICunha & Smltl] i2Q06ll . Stellar pai'ameters from iRich & OrigU^ i20Q5h . 
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(Australia). More inf ormation on the instrument can be found in 
iRodgers et al.l ( Il988h . Exposure times ranged between 200 and 
1200 s. The precise wavelength range of the optical spectra varies, 
but the part between 666 - 1000 nm is always covered. This range 
includes prominent TiO and VO bands that are very sensi tive to 
the effective temperature of the star (see e.g. lReinersll2005h . Fur¬ 
thermore, the general slope of the spectra can constrain the star’s 
temperature. 

For the Spitze r sample stars identified as non-variable by 
IVanhollebekd | |2007L Sect. lS.l.Tt . low-resolution spectra that were 
obtained several years before the CRIRES run were used. These 
older spectra have a dispersion of 0.41 nm/pixel. The sample stars 
identified as variable were observed with the DBS between 16 and 
18 days after the CRIRES run in which the Spitzer sub-sample 
was observed. These spectra have a somewhat higher resolution of 
0.187 nm/pixel. We may safely assume that the temperature derived 
from the optical spectra is representative for the temperature at the 
time of the CRIRES observations: stars with short pulsation peri¬ 
ods and low amplitudes do not exhibit a large temperature change 
over a pulsation cycle, while for the stars with longer periods and 
larger amplitudes the time delay between CRIRES and DBS ob¬ 
servations is small compared to the pulsation period. For the Plant 
stars, the CRIRES and DBS observations were done at indepen¬ 
dent phases, hence here the temperature determined from the low- 
resolution spectra may deviate somewhat more from that at the ac¬ 
tual time of CRIRES observations. All sample stars were identified 
to be oxygen-rich from the low-resolution spectra, i.e. they all ex¬ 
hibit the characteristic TiO bands. Also the Tc-rich stars have very 
prominent TiO bands, the ZrO bands are present but not dominant. 
In the observed spectral range, C-stars would be identified by in¬ 
tense bands of the CN molecule. 

Before temperature determination, the optical low-resolution 
spectra were correct ed for interste ll ar red dening. The same ex¬ 
tinction values as in [Wnhollebeta j 20()7) w e re adopted, w hich 
stem from the maps of Tschultheis et alj ( Il999h . ISum 1 ( |2004 . and 
iMarshall et Zl l l2006h . These values were applied to establish the 
extinction as a fun ctio n of wavelength, usin g the relations given in 
lO’ Donnell jl994h and ICardelli et alj h989l) . The spectra were di¬ 
vided by this extinction curve to obtain the de-reddened spectrum. 
Test runs show that the resulting temperatures are quite insensitive 
to the extinction value. 

The main method to derive temperatures from the optical spec¬ 
tra involves fitting of series of model spectra to the observed spec¬ 
tra. A series of COMARCS model atmospheres and spectra (see 
Sect. 13. 2t was calculated, with temperatures between 2600 and 
3800 K in steps of 100 K, loggr[cms“^] = 0.0, 1 Mq, microtur¬ 
bulence 3 km s“', solar C/0 ratio, and two metallicities: 1.0 Zq and 
0.3 Zq. Extensive lists of atomic and molecular lines, in particular 
of TiO and VO, were taken into account in the calculation of the 
spectra. In addition , we used a grid of PHOENIX model spectra 
jHusser et al. I [20l3h with temperatures between 2000 and 3900 K 
in steps of 100 K, solar metallicity, and otherwise identical param¬ 
eters as the COMARCS grid. We fitted the model spectra to the 
observed spectra with a modified method in the form 


_ 1 ^ - /model., ([Fe/H]))" 

k (/obV + /model,, ([Fe/H])) 


The wavelength range between 700 and 930 nm was used to find the 
temperature at which the model spectra best fit an observed spec¬ 
trum. The lower wavelength limit was chosen to avoid the noise- 
dominated blue part of the observed spectrum, whereas the upper 


limit was chosen to avoid a molecular band of ZrO that might be 
present in the stars that had experienced 3DUP but not in stars that 
had not. This band could potentially have a systematic impact on 
the derived temperature. For the non-variable stars in th e sample, 
the (V - K)q vs. Teff calibration of ivan Belle et al.l ( Il999h was used 
to derive another temperature estimate. This photometric tempera¬ 
ture estimate agrees well with that derived from fitting the optical 
spectra. The final adopted temperature listed in column 6 of Table[T] 
was chosen based on a visual inspection of the observed spectra and 
of the model fits to them. We also tried to derive tem peratures from 
the spectral parameters defined in lFluks et al.l i ll994 their Eq. 17 - 
19), but found that these estimates sometimes vary considerable for 
a given star. The narrow bandpasses lie at the bottom of molecular 
bands that are formed in the very outer layers of the star, which are 
most affected by dynamic effects (pulsation, mass loss). 

The determined temperature may depend on the adopted 
metallicity of the model spectra. Conversely, the metallicity that we 
derive from the high-resolution spectra depends on the model input 
temperature. An independent determination of effective tempera¬ 
ture and metallicity without any assumptions is very difficult for 
our sample stars. In the above procedure we find relatively mod¬ 
est differences between the temperatures derived using the models 
with 1.0 Zq and 0.3 Zq, respectively. The impact of the assumption 
on the metallicity for temperature determination seems limited and 
we are thus confident that we can derive reliable abundances from 
the high-resolution spectra. In any case, the CO lines are quite in¬ 
sensitive to the effective temperature. 


3.1.3 Luminosity 


The luminosity of each sample star was determined by inte¬ 
grating its spectral energy distribution (SED). To construct the 
SEP, we collected I- band photometry from the DENIS catalogue 
( lEpchtein etalj |l997[). cycle-a verag ed J, H, K, L^ h band photom¬ 
etry from Wood et ^ ( I l998h and Vanhollebeke |2007|). Soitzer 
IRAC and MIPS 24 photometry dUttenthaler et alj201(]l : lHinz et al.l 
I2OO9I) . as well as photometry f rom the WISE jW right et a l]|201(]h . 
ISOGAL I Omoniet^ l20()3t) . Akari dlshihara et al. 2010h . and 
MSX l lEgan et al.ll200 ^ surveys. Because for most stars the bulk 
of the stellar flux is emitted between the J and L bands and be¬ 
cause cycle-averaged fluxes in these bands are used, we are con¬ 
fident to have very precise luminosities at hand for our sample 
stars. The photometry was de-reddened for interstellar extinction 
(see Sect. 13.1.^ and converted to flux values (Jansky). The whole 
SED was integrated numerically by connecting the data points by 
straight lines, and an extrapolation to the origin (y, /) = (0,0) was 
included. To convert this apparent flux to a luminosity in solar units, 
a distance modulus of 14"’5, corresponding to a distance of 8.0 kpc 
to the bulge, was assumed. A typical example of a de-reddened 
SED is shown in Fig.[T] 

For stars with a measured pulsation period, we corrected these 
luminosities by forci ng them to fall on one of the period - K- 
magnitude relations of llta et aT] ( l2004h . in an attempt to correct for 
the distance scatter within the bulge (red points in Fig.[3j. A differ¬ 
ence in the distance modulus between the Large Magellanic Cloud 
and the bulge of 4?0 was assumed. A log P — Kq diagram (Fig. O 
reveals that the long period stars probably follow sequence C (fun¬ 
damental mode or Mira-like pulsators), whereas the shorter period 
stars probably follow sequence C’ (first overtone or semi-regular 
pulsators). We assume that the transition occurs at a period of 160 d 
(log(P) ~ 2.2). Both the uncorrected and corrected luminosities are 
listed in columns 7 and 8 of Table[T] 
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Figure 1. Typical example of a de-reddened SED of a sample star 
(J174157.6-282237). The source of the various data points is indicated in 
the legend. The luminosity of this star was found to be 3342 Lq. 



Figure 2. log P — Kq diagram of our sample stars. Black symbols are the 
Spitzer sample stai's, filled red diamonds the Tc-rich Plaut stars, open red 
diamonds the Tc-poor Plaut stars. The solid lines are the relation s C (fun¬ 
damen tal mode pulsators) and C’ (first overtone pulsators) from llta et alJ 

<2004 . 


This correction would lead to very high luminosities (8000 Lq 
and more) and subsequently unrealistically high current masses 
(3-5Mq, see next Section) for the five Spitzer sample stars with the 
longest pulsation periods, four of them being the LAVs. Because 
their very red {J — K)q colours (up to 2.29) and considerable mid-IR 
excesses suggest that their K-band magnitude s are diminished by 
strong circumstellar extinction (cf. Fig. 20 in IVassiliadis & Wood 
Il993h we did not apply the distance correction to their luminosi¬ 
ties. The long-period Plaut stars do not show such extreme colours, 
which is why the correction was applied nevertheless. 

The Plaut stars preferentially lie above sequence C in Fig. 
whereas the Spitzer stars seem to lie preferentially below the 
sequences. The location of the Plaut stars can be understood 
as a result of selection effects that favoured stars on the near 
side of the bulge. They were discovered by photographic plate 
photometry dPlau 3ll97lh and selected for optical spectroscopy 
(lUttenthaler et al. 20071) . The location of the Spitzer stars is more 
surprising, at the moment we do not have a convincing explanation 



Figure 3. Luminosity vs. effective temperature diagram of the sample stars. 
Stars whose luminosities were not corrected for the bulge depth scatter due 
to the lack of a known pulsation period are plotted as black symbols, those 
whose luminosities were corrected are plotted as red symbols. The Plaut 
stars are plotte d with the same sym bols as in Fig. |2] The blue graph is an 
isochrone from iBressan et al]j2012h with 5.0 Gyr age, solar metallicity, and 
including dust formation on the AGB. 


why they are preferentially at distances farther than the Galactic 
centre. 

Figure [2 shows a Hertzsprung-Russell diagram with inverted 
temperature axis of the sa mple stars. The solid blue line in that 
figure is an isochrone from iBressan et all ( 1201 2h with 5.0 Gyr age, 
Z = 0.017 (solar metallicity), and including dust formation on the 
AGB. Stars on the AGB of this isochrone have a current mass of 
~ 1.29Mq. The RGB tip of this isochrone is at ~ 25OOL0. Most of 
the sample stars define a clear trend of increasing luminosity with 
decreasing temperature, as may be expected from stars distributed 
along the AGB. Because temperature and luminosity of the stars 
were derived from completely independent data, it can be excluded 
that a systematic effect produces such a trend. The isochrone covers 
well the distribution of of our sample stars and reproduces the trend 
of increasing luminosity with decreasing temperature. The Tc-rich 
Plaut stars are shifted away from the isochrone to somewhat higher 
temperature and/or higher luminosity. At lower metallicity, stars 
on the AGB have higher temperatures at a given luminosity. As 
will be shown later, the two Tc-rich stars for which we are able 
to determine the metallicity are relatively metal-poor. This could 
explain their displacement from the location of the other stars. 


3.1.4 Current mass and surface gravity 

Current pulsation masses for the variable stars with a determined 
period were computed. The linear pulsation models used here have 
been applied to LMC cluster variables, where they yield masses 
that are in good agreement with those expected fr om the cluster 
ages determined by main-sequence turn-off fitting jKamath et al.l 
iTOlOh . Li near pulsation periods wer e computed with the code de¬ 
scribed in iLebzelter & Wo^ (120071) . except that the low temper- 
ature opacities were gen erated from the web site described in 
iMarigo & Aringeil ( l2009h . Metal and helium mass fractions of 
Z = 0.016 and Y = 0.28 were used. A grid of AGB models with 
600 < L/Lq < 10000 and 0.7 < MJMq < 5.0 was made and the 
linear fundamental and first overtone periods calculated. Analytic 
fits to the models were made in the form 
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log ^ = 0.324 - 1.699p + 3.050r + 1.007p^ + 1.891r- 
and 


•2.515pr 


■ 0.894r^ + 1.524pr 


log — = -0.143 - 2.148p + 3.271r - 0.687^^ 

Mq 

for the fundamental and first overtone modes, respectively, where 
p = log P-2 and r = 0.5 (log L/Lq - 4 log Teff) + 5.3. These fits to 
the pulsation mass are accurate to better than 7% for all computed 
models. Stars with L/L© < 2000 could be RGB stars rather than 
AGB stars. Similar fits were made to RGB models with luminosi¬ 
ties 600 < L/T© < 2000 and 0.6 < T//M© < 1 . 0 , the difference 
between the AGB and RGB models being the core mass. The core 
masses for the AGB models w ere assumed to be given by the for¬ 
mula in IWood & Zarrol (Il98lh while RGB masses were obtained 
from the stellar models of Girardi et al.l l l2000h . The RGB fits are 


log-^ =0.351 
Mq 

and 


1.900p + 3.397r + 1.303p^ + 2.917r- - 3.750pr 


log-= -0.165 - 2.521p + 3.682r - 1.256p^ - 0.952r^ + 2.343pr 

Mq 

for the fundamental and first overtone modes, respectively, with an 
accuracy better than 2 %. 

Using the values of P, LILq (corrected with the help of pe¬ 
riod - K-magnitude relations), and Teft given in Table [T] pulsation 
masses were computed from the above formulae. It was assumed 
that stars with P < 160 d are first overtone pulsators and those 
with P > 160d are fundamental mode pulsators (cf. Fig. |2j. The 
pulsation masses are listed in column 9 of Table [T] Since for stars 
with L/Lq < 2000 the difference between RGB and AGB pulsa¬ 
tion masses is small, only the AGB pulsation mass was adopted. 
The masses of the five Spitzer sample stars with the longest pulsa¬ 
tion periods have to be considered the least certain. 

With these current masses estimated from pulsation theory, the 
surface gravity of the sample stars was calculated using the “clas¬ 
sical” relation: 


-log(^) + log 

adopting Tq = 5770 K and loggr© = 4.44. A few stars turned 
out to have unrealistically low masses (M < O.IMq), i.e. masses 
less than the carbon-oxygen core plus a thin envelope that is ex¬ 
pected to be left when the star leaves the tip of the AGB. To obtain 
more realistic log o' values for them, a current mass of O.IMq was 
adopted. For stars without a known pulsation period, a current mass 
of M = 0.85M© was adopted as a first guess. This is a reasonable 
assumption for a supposedly old population. However, for most of 
these stars, this lead to a very low surface gravity, lower than for 
what the COMARCS atmospheric models converged (Sect. Ira . 
These stars are also outliers in a log g vs. luminosity diagram. Thus, 
for them we used the most extended model atmospheres that were 
available (i.e. lowest log<;). These adopted log 9 values are well 
within the range expected from the general trend of log g vs. lu¬ 
minosity formed by the other stars. Column 10 of Table[T]lists the 
surface gravities log g finally adopted for the model atmospheres. 


logg* =logg© + 41og(^ 



transfer routines from lNordlundl i ll9841) and assumes local thermo¬ 


dynamic equilibrium (LTE). The opacities are taken from tables 


prepared by the COMA code jArinaei 

|2000t). Recent COMARCS 

models are described injAringer et al. i 

2005i. 


3.3 Measuring radial velocities and abundances 

Radial velocities (RVs) of the stars were determined from the nu¬ 
merous CO Av = 2 lines in the K-band setting in grating order 
24 by cross-correlation with a generic synthetic model spectrum. 
These lines are quite insensitive to abundances and stellar parame¬ 
ters, hence the choice of the model has a negligible impact on the 
measured radial velocities. The typical uncertainty on the RV is 
~ 1 km s“', but may be better for the hotter, less variable sample 
stars. The radial velocities are summarised in column 2 of Table|2] 
COMARCS model atmospheres adopting the main stellar pa¬ 
rameters (7’eff,logg) derived in Sect. l3.1l were calculated for each 
star, spanning a grid of various general abundance combinations. 
In a first try, we established a three-dimensional grid in [Fe/H], 
C/0, and [a/Fe], where the a-elements are O, Ne, Mg, Si, S, Ar, 
Ca, and Ti. [Fe/H] varied between -1.5 and +0.5 in steps of 0.5, 
for C/O the values 0.15, 0.30, 0.60, and 0.90 were adopted, and 
[o'/Fe] varied betwe en -0.20 and +0.40 in steps of 0.20. The solar 
abundances listed in ICaffau et all ( l2008h were adopted as reference 
scale. The solar metallicity on that scale is Z© = 0.0156 and the 
solar C/O ratio is 0.55. A generic microturbulence of 2.5kms“', 
a value t ypical for evolved giants, was adopted for the model cal¬ 
culation l lCunha et al.ll2007l ]. Based on these model atmospheres. 


synthetic spectra were generated for the wavelength range of the 
observed H-band spectra. Important line lists that wer e taken into 


account were the atomic line list of R^deet^ 

2 OIOI 1 . the CN as 

well as the H 9 O line lists fromlldrsenseiTi 19^ 

, the OH list from 

the HITRAN 2004 comnilation 1 Rothman et alj 

200.5b. and the CO 

list fromiGoorvitch & Chackerian| 1 199411. For both the CN and CO 


line lists the wavelengths of so me transitions were corrected with 
the help of line identifications of lHinkle et al.l ( ll995t) . The synthetic 
spectra with an original resolving power of R = 300 000 were con¬ 
volved with a Gaussian profile to that of the observed spectra at 
R = 50000. We developed an IDL routine that interpolated be¬ 
tween these grid spectra and fitted the interpolated model spec¬ 
trum to the observed one of each star. This interpolation was done 
for every point of the wavelength vector. The standard IDL rou¬ 
tine amoeba. pro was employed to find the interpolated abundance 
combination that achieved the best fit to the observed spectrum, by 
minimising a^^ criterion. 

The initial aim was to derive the abundance of oxygen, the 
most abundant a-element, from lines of the OH molecule that are 
prominently present in the H-band. However, it turned out that the 
OH lines present in the spectra are not sensitive enough to use them 
as reliable abundance indicator of O. In a consistent spectral syn¬ 
thesis (i.e. the same chemical abundances are adopted for the model 
atmosphere and spectral synthesis calculation), the strength of the 
OH lines is almost unchanged, even in a wide range of O abun¬ 
dances. Such a low sensitivity of molecular line strength to abun¬ 
dance is known to occur in cool stellar atmospheres and can be 
explained by the feedback of the abundance of key elements, most 


3.2 Model atmospheres 

Model atmospheres for our sample stars were established with the 
help of the COMARCS code. COMARCS is based on a revised 
version of MARCS JjOrgensen et alJl992l) with spherical radiative 


^ This value is slightly lower than the microturbulence used in the grid used 
for temperature determination, but the difference has virtually no impact on 
the results. 

^ This H 2 O line list was chosen due to its completeness. 
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Table 2. Radial velocities and abundances measured from the CRIRES spectra. Meaning of the columns: object 
name; heliocentric radial velocity; overall [Fe/H] metallicity; C/0 ratio; carbon isotopic ratio; relative abun¬ 
dances of Al, Si, Ti, and Y. Abundances of stars for which no acceptable fit was achieved are omitted (...), 
abundances followed by a colon are uncertain (poor fit). 


Object name‘s 

I^^helio 
(km s-^ 

[Fe/H] 

C/0 

‘2C/‘3C 

[Al/Fe] 

[Si/Fe] 

[Ti/Fe] 

[Y/Fe] 

BMB 78“ 

-59.7 

-0.47 

0.46 

9.9 





BMB 289“ 

-65.2 

-0.05 

0.57 

4.0: 





BMB 289'’ 


-0.21 

0.67 

4.6: 





BD-012971'’ 

-119.6 

-0.68 

0.30 

9.5 

+0.40 

+0.21 

+0.60 

-0.07 

J174117.5-282957 

159.4 



9.3: 





J174123.6-282723 

66.5 

-0.07 

0.39 

18.2 

+0.05 

-0.04 

+0.36 

-0.02 

J174127.3-282851 

58.7 








J174127.9-282816 

-166.0 

0.22 

0.22 

14.5 

-0.03 

-0.10 

+0.17 

+0.11 

J174128.5-282733 

15.0 

-0.21 

0.25 

14.6 





J174139.5-282428 

102.3 



4.8: 





J174140.0-282521 

-89.7 



11.3: 





J174155.3-281638 

177.8 



11.9: 





J174157.6-282237 

-101.7 

-0.18 

0.42 

5.5 

+0.09 

+0.16 

+0.50 

+0.19 

J174203.7-281729 

-234.2 








J174206.9-281832 

56.2 








J174917.0-293502 

18.1 

+0.01 

0.50 

19.4 

+0.18 

-0.07 

+0.20 

0.00 

J174924.1-293522 

2.2 

-0.10 

0.23 

18.3 





J174943.7-292154 

-9.3 

-0.16 

0.44 

22.2 

-0.05 

-0.02 

+0.23 

-0.18 

J174951.7-292108 

-72.4 

-0.06 

0.51 

17.2 





J175432.0-295326 

15.5 



17.7: 





J175456.8-294157 

107.1 



17.4: 





J175459.0-294701 

198.6 



9.8: 





J175515.4-294122 

171.0 

+0.01 

0.44 

16.5 





J175517.0-294131 

-55.0 

-0.08 

0.34 

16.9 





J180238.8-295954 

103.2 

-0.17 

0.45 

16.0 

+0.25 

+0.23 

+0.43 

+0.13 

J180248.9-295430 

143.3 

-0.15 

0.48 

20.7 

+0.21 

+0.01 

+0.41 

+0.19 

J180249.5-295853 

^2.3 

-0.30 

0.29 

12.8 





J180259.6-300254 

^0.3 

-0.33 

0.26 

22.3 





J180301.6-300001 

58.4 

-0.05 

0.46 

19.3 





J180304.8-295258 

33.7 

0.07 

0.33 

9.9 





J180305.3-295515 

-21.4 

-0.27 

0.30 

9.8 

+0.17 

+0.18 

+0.44 

+0.15 

J180305.4-295527 

27.5 

-0.18 

0.26 

15.5 

+0.16 

+0.18 

+0.39 

+0.02 

J180308.2-295747 

14.9 



11.0: 





J180308.7-295220 

13.7 

-0.12 

0.51 

27.9 

+0.23 

+0.11 

+0.38 

+0.02 

J180311.5-295747 

44.8 








J180313.9-295621 

-96.4 

-0.26 

0.29 

11.0 

+0.15 

+0.04 

+0.23 

-0.11 

J180316.1-295538 

62.5 

-0.15 

0.47 

18.4 





J180323.9-295410 

^2.4 

-0.21 

0.37 

18.8 

+0.11 

+0.01 

+0.21 

-0.07 

J180328.4-295545 

-113.7 

-0.17 

0.32 

7.5 





J180333.3-295911 

6.6 

-0.56 

0.39 

9.7 





J180334.1-295958 

24.0 

-0.11 

0.47 

16.1 

+0.31 

+0.27 

+0.35 

+0.09 

Plaut 3-45 

16.3 








Plaut 3-100 

-22.9 








Plaut3-315 

-61.4 








Plaut 3-626 

-67.2 

-0.85 

0.71 






Plaut 3-794 

-50.1 








Plaut 3-942 

-80.1 

-0.56 

0.75 

19.6: 





Plaut3-1147 

-2.3 








Plaut 3-1347 

63.0 









^ Abundances ado pting the stellar param eters from ICunha & SmithI bOOfil) . Abundances adopting the 
stellar parameters from lRich & Origlia i2005l) . 
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Figure 4. Synthetic spectra based on COMARCS atmospheres with Feu = 
3500 K, log 3 = +0.25, [Fe/H] = +0.16, and [ff/Fe] ratios as indicated in 
the legend. The model atmospheres and synthetic spectra were consistently 
calculated with these abundances while keeping the carbon abundance con¬ 
stant. 


importantly O, on the atmospheric structure due to their impact on 
molecular opacities. The atmospheric structure may change just in 
such a way as to leave OH line strengths unchanged. This point is 
illustrated in Fig. |4l which shows a zoom-in on synthetic spectra 
around three OH lines just hlue-ward of the '^CO 3-0 hand head. 
The CO and CN lines become weaker with increasing oxygen abun¬ 
dance because the chemical equilibrium is changed to the disfavour 
of these molecules, whereas the OH lines remain unchanged with 
O abundance. This conclusion from our models was also checked 
and c onfirmed with models from the PHOENIX code dHusser et al.l 
l2013h . Hence, from the OH lines available in our spectra the O 
abundance cannot be derived with a sufficient precision in cool, 
M-type giants, and all previous such determinations that used an in¬ 
consistent calculation have to be taken with caution. Only for hotter 
stars such as Arcturus (T^s > 4000 K) the OH lines are sensitive to 
the O abundance. 

Because it is not possible to constrain the a-abundance from 
the H-band spectra, the original 3D grid of model atmospheres had 
to be reduced to a 2D grid. Here we make the crucial assumption 
that [or/Fe] is a function of [Fe/H] as follows: [a/Fe] = +0.40 
at [Fe/H]=-1.5, -1.0, and -0.5, followed by a linear decrease 
of 0.4dex/dex up to [Fe/H]=+0.5. T his run closely appr oximates 
that found from dwarf and sub-giant dBensbv et alJboO. see th eir 
Fig. 25) as well as giant bulge stars (e.g. Alves-Brito et al.f2010l) . 

In addition to fitting the two detector chips from the H-band 
setting, chip 2 from the additional K-hand setting grating order 26 
was included when available. This chip was chosen because it con¬ 
tains a number of metal lines that are not too strong, but only a rel¬ 
atively small number of CN lines. The metal lines should provide 
for a better constraint on the metallicity of the stars, while the CN 
lines are also sensitive to the C/0 ratio. The change in [Fe/H] and 
C/O implied by including this additional piece of K-hand spectrum 
is relatively small. For the 13 science targets with reliable abun¬ 
dance determination (i.e. excluding J174128.5-282733, J180308.2- 
295747, and J180311.5-295747), the average change is 0.00 dex, 
with a standard deviation of 0.07 dex. 

A map ofx^ in the ([Fe/H], C/O) plane shows that the^‘ val¬ 
ues form a relatively flat valley that runs diagonally through that 
plane: A similarly good fit as at the minimum is achieved at a lower 


(higher) metallicity, if at the same time the C/O ratio is increased 
(decreased). By including the additional chip from the K-band, the 
walls of the “valley” become steeper along the [Fe/H] axis, which 
means that the metallicity is better constrained for the stars with 
these additional observations. 

A first test of the analysis methods was done with synthetic 
observations, i.e. synthetic spectra that were created with model 
atmospheres that have abundance combinations that are between 
the model grid. Gaussian noise was added to the spectra to simulate 
an S/N ~ 40. Also, a slightly altered line list was used for the 
synthetic observations to simulate absorption lines missing in the 
models. It turned out that the fitting routine tended to converge to 
a metallicity that was slightly (~ 0.10 dex) lower than the input 
metallicity, and a C/O ratio that was slightly higher than the input 
C/O. The results improved considerably when wavelength points 
were neglected from the fitting procedure whose flux deviated too 
much (~ 0.10 of continuum flux) from the observed flux at a given 
iteration step. Hence, this clipping technique was also adopted for 
the analysis of the real spectra. 

The results of the fitting procedure were checked by visually 
inspecting a number of metal lines that were found to be sensitive 
to [Fe/H], but not to C/O (i.e., free of CN and CO blends). This 
included the following lines: Fe 1549.46, Fe 1553.60, Fe 1553.85, 
Fe 1557.10, and Ni 1561.00 nm. If these lines were not well repro¬ 
duced by the automatic fitting procedure, the fitting was repeated 
manually to reproduce the strength of these lines. Subsequently, 
the fitting procedure was run again to measure the C/O ratio of the 
star with the metallicity fixed to the value derived manually. 

A typical example of a spectral fit is shown in Fig.|5] No dy¬ 
namic effects such as line splitting, emission features, etc., are seen 
in the warmer (Teff > 3100 K) stars, even in the variable ones. Most 
of the large residuals (observed minus calculated flux) do not vary 
in wavelength from star to star; t hey ar e most probably caused by 
missing line data, cf. iRvde et alj l l2009l) for some unidentified fea¬ 
tures. 

We may thus assume that the ability to model the spectra of 
our sample stars is mainly limited by the (in)completeness of avail¬ 
able line lists, except for the coolest stars (T^g < 3000 K), for which 
the hydrostatic models are not realistic and water starts to play 
an important role. At this low a temperature, the metal lines lose 
their power to constrain the metallicity and the results for these 
cool stars are not reliable. Thus, abundances for stars cooler than 
Teff = 3000 K are not reported in Table[^ 

Most notably, it is not possible to model the spectra of the 
four LAVs in the Spitzer sample, even though they were observed 
close to minimum light when their spectra are supposed to be least 
complex. A m o del sp ectrum using only the line list of water by 
iRothman et alj l l2010n was calculated, based on one of the cool 
models (T^g = 2800 K) to compare it to observed cool star spectra. 
Many of the strong lines in the observed spectra have a counterpart 
in the synthetic water spectrum, suggesting that water is indeed a 
strong contributor to the spectrum of these stars. This comparison 
also shows that the match in wavelength between observed spec¬ 
trum and water line list is good at this high a resolution. 

The abundance determination was however successful for the 
two warmest Plaut stars, even though they are large amplitude vari¬ 
ables. 

We note that our analysis represents a significant effort. The 
computation of a single opacity table with one abundance combina¬ 
tion takes 19 hours on a fast, modem desktop PC. The construction 
of the model atmospheres requires significant human interaction 
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Figure 5. Typical example of a spectral fit achieved in this work, for star J180313.9-295621. The upper, middle, and lower panel show chip 2 and 3 of the 
H-band setting (order 36), and chip 2 from the K-band setting with order 26, respectively. The upper sub-panel in each panel shows in black the observed 
spectrum and in red the best-fitting interpolated synthetic spectrum. The lower sub-panel in each panel shows the residuals (observed - calculated flux) of the 
fit. This star was found to have Teff = 3400 K, log g = +0.05, [Fe/H] = —0.26, and C/O=0.29. 
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Figure 6. Fit to the *'^CO 2-0 R76 line of the star J174123.6-282723. 
The solid red line shows the synthetic spectrum interpolated to the best-fit 
'^C/'^C=18.2. The dotted lines show the synthetic spectra calculated with 
'^C/'^C=50, 12, and 3.5 (from top to bottom). 


because the models do not easily converge at low logg. A more 
detailed analysis is beyond the resources available to us currently. 


3.4 The carbon isotopic ratio and metal abundances 

To determine abundances of metals as well as the carbon isotopic 
ratio, we established for each star a tailored model atmosphere 
adopting the stellar parameters derived in Sect. 13.11 and assuming 
the general abundance pattern ([Fe/H], C/O) as found from the H- 
band spectra (and chip 2 of order 26, where available). 

The carbon isotopic ratio was derived from lines of *'^CO lo¬ 
cated in chip 1 of grating order 24. An inspection of the optical 
depth of line cores as a function of gas temperature reveals that 
the low-J (R19, R22, R24, R26) *^CO 2-0 lines are formed rela¬ 
tively far out in the atmosphere where dynamic effects may cause 
significant deviations of the model atmosphere from the real atmo¬ 
sphere. The analysis was thus based on the high-J 2-0 lines (R76, 
R78, R80, R81). These lines are also more sensitive to changes in 
than the low-excitation lines. Spectra with ratios 

of 3.5 (equilibrium value of the carbon isotopic ratio in the CN cy¬ 
cling), 7, 12, 20, 30, 50, and 90 (solar value) were synthesised based 
on the tailored model atmospheres. The isotopic ratio of the stars 
cooler than 3000 K, for which no general abundance pattern could 
be determined, was measured with a generic abundance combina¬ 
tion of [Fe/FI]=0.0, [ff/Fe] = -(-0.20, and C/O=0.30. Their isotopic 
ratios are more uncertain and are marked by a colon in Table [2] 
Again, an interpolation and fitting routine was employed to deter¬ 
mine the value of isotopic ratio that best fits the observed spectrum 
in the mentioned lines. Figure[^shows a typical example fit to the 
‘^CO 2-0 R76 line of the star J174123.6-282723. Note the good fit 
also to the *^CO 4-2 R64 line at the right hand side of the figure, 
which gives additional confidence in the C/O ratio that was derived 
for this star (0.39), though the CO Av = 2 lines in the K-band are 
generally less sensitive to C/O than the CO Ay = 3 lines in the 
Fl-band. 

We also aimed at deriving abundances of individual metals us¬ 
ing equivalent width (EW) analysis. Weak enough (log(EWobs/'l) < 
-4.5) but unblended lines of Al, Si, Ti, and Y were found in the 
spectra in grating orders 26 and 27. Y is a heavy element produced 
in the slow neutron-capture (s-) process, an increase in its abun- 


Table 3. Line list used in the determination of the metal abundances. Mean¬ 
ing of the columns: chemical symbol of the element; wavenumber of the 
transition; oscillator strength gf; excitation energy of the lower level in 
10‘*cm-‘. 


Element 

O’ 

(cm-‘) 

gf 

(10‘*cm-‘) 

All 

4713.8730 

3.291e-l 

4.13194 

AH 

4713.8810 

3.291e-l 

4.13194 

Si I 

4779.4609 

3.104e-b0 

5.42569 

Til 

4589.4991 

1.256e-l 

1.41066 

YI 

4702.2860 

l.OOOe-1 

1.15337 


dance can indicate recent or ongoing 3DUP activity in an AGB star. 
Astrophysical g f values of the lines were established with the Arc- 
turus spectrum dlTinkle e t alJ 199 5|). adopting the abundan ces for 
Al, Si, and Ti reported in Ramirez & Allende Prietol (1201 ih . and a 
scaled solar abundance f or Y. Other atomic parameters were taken 
from the VALD database dKunka et al.l200(i) . The resulting line list 
is reported in Table [3] Synthetic spectra with varying abundances 
of the individual elements were calculated from the tailored model 
atmospheres, from which the EWs of the respective lines were mea¬ 
sured. The abundance in each star was then calculated by interpo¬ 
lating in the log(EWsy„th/d) vs. loge(X) plane to the observed re¬ 
duced EW (log(EWobJ4)). 


3.5 Abundance error estimate 

To estimate the formal errors on the measured abundances, uncer¬ 
tainties in four quantities were taken into account: temperature, sur¬ 
face gravity, continuum placement, and microturbulence. 

An estimate of the temperature uncertainty may be derived 
from the scatter that most stars exhibit around the general trend 
of decreasing temperature with increasing luminosity in the HRD 
(Fig.|3). If the LAVs and the obviously deviating Plant stars are 
excluded, the remaining sample stars follow a linear least squares 
fit in Teft vs. uncorrected luminosity with a standard deviation of 
181 K. For the stars with known pulsation period the scatter re¬ 
duces from 175 K before correction to 119 K after correction. Note 
that the relation between temperature and luminosity may have an 
intrinsic scatter and curvature (as indicated by the isochrone in 
Fig.H, and that the temperature is not pulsation-cycle averaged. 
We may thus assume that the uncertainty of the temperature deter¬ 
mination is not larger than this scatter. An uncertainty of ±100K 
was adopted. 

To estimate the uncertainty on the surface gravity (log g), we 
adopted the temperature uncertainty, an uncertainty on the luminos¬ 
ity of 10%, and on the pulsation period of 1%. These uncertainties 
were applied to the formulae given in Sect. l3.1.4l From this, a typ¬ 
ical uncertainty of A log g = ±0.12 was derived. 

The continuum placement was changed by ±1% and the mi¬ 
croturbulence ^ by ±lkms“*. The latter seemed a reasonable un¬ 
certainty from the CO lines in the K-band, which are quite sensi¬ 
tive to To determine the abundance uncertainty, the model atmo¬ 
spheres and synthetic spectra were recalculated, varying the param¬ 
eters one at a time by the given uncertainties, and the spectral fitting 
was repeated for a few representative stars. The total uncertainty on 
the abundances was calculated by combining the individual sources 
of error in quadrature. 

The general abundance pattern ([Fe/H], C/O) turned out to be 
quite robust. The metallicity is most sensitive to temperature and 
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whereas the C/0 ratio is most sensitive to changes in log g and ^ 
(the CO lines in the H-band are quite sensitive to log g). Changes 
in [Fe/H] due to temperature uncertainty vary between 0.01 and as 
much as 0.21 dex, but usually they are ~ 0.06 dex. Typical com¬ 
bined uncertainties on [Fe/H] and C/0 are 0.10 dex and 0.08, but 
C/0 is more robust at high values (±0.05 at C/O=0.50). The C/O ra¬ 
tio was found to change by about 0.1 if log g is varied by ±0.3 dex. 
This means that changing the stellar mass, and hence log g, within 
any reasonable limits does not significantly change our results and 
conclusions on the C/O ratio of the stars. 

The parameter having the largest impact on the *^C/*^C ratio 
is the C/O ratio. Changing it by 0.10 can change the isotopic ratio 
by as much as ~ 30%. The errors are correlated: An overestimate 
(underestimate) of C/O will also lead to an overestimate (under¬ 
estimate) of *-C/*'^C, because with initially stronger (weaker) CO 
lines it takes less (more) *^C to get the same line strength. The sur¬ 
face gravity and continuum placement have some influence on the 
'^C/'^C ratio, but is limited to ~ 10% and ~ 7%, respectively. We 
estimate a combined uncertainty of ± ~ 33% on the '^C/'^C ratio, 
comparable to literature values. 

Typical uncertainties for the elemental abundances are ±0.11 
for [Al/Fe], ±0.25 for [Si/Fe], ±0.13 for [Ti/Fe], and ±0.16 for 
[Y/Fe]. 

3.6 Comparison stars 

To test our tools for abundance determination, we applied them to a 
few comparison stars. The goal of this exercise was not to check lit¬ 
erature results, but rather to compare abundances derived with our 
tools to those from other studies and to identify possible systematic 
differences. 

The most important comparison object in studies of cool 
stars is of course Arcturus. A grid of model atmospheres span¬ 
ning the same range of [Fe/H] and C/O values as for the science 
tar gets was ca l culate d, adopting the stellar parameters determined 
bv iRvde et alj (|201(]|) . which are: = 4280 K, log^ = 1.7, and 

<fmicro = 1.7kms“*. Synthetic spectra were generated and fitted to 
the ob served, high-resolution spectrum provided by iHinkle et al.l 
( Il995h in the same spectral range in which the science targets were 
observed with CRIRES. The fitting routine converged at a metal- 
licity of [Fe/H] = -0.50, which a grees very well wit h the value of 
[Fe/H] = -0.53 ± 0.05 derived bv iRvde et al] (l201(ll) . Our C/O ra¬ 
tio of 0. 32 is somewhat higher than 0.24 as found by these authors, 
whereas [Rich & Origlial (l2005h determined C/O=0.28 in Arcturus. 
Since Arcturus is a K-type giant whose atmospheric structure is 
less influenced by molecular opacity than that of M-type giants, 
we also performed the 3D fit introduced in Sect. l3.3l With this ap¬ 
proach the following abundance pattern for Arcturus was found: 
[Fe/H] = -0.51, C/O=0.37, and [n/Fe] = ±0.31. The a-abundance 
ratio is in e xcelle nt agreement with the value of ±0.33 found by 
iRvde et ^ ( 1201 Oll from an average of the abundance of O, Si, S, 
and Ti. We also modelled the observed CRIRES spectra of the 
comparison stars BMB 78, BMB 289 (both located in the Galactic 
bulge), and BD-012971 (a solar neighbourhood star). The stellar 
parameters of these stars were not deter mined as for the science tar- 
gets, rather the parame ters provided bv ICunha & Smitl] ( l2006l) and 
iRich & Origlial ( l2005h . respectively, were adopted (see Table [T]l. 
The results of our abundance determination as well as the literature 
values are summarised in Table|4] 

The differences to the literature values vary from star to star. 
Our [Fe/H] is lower than the literature [Fe/H] value by more than 
0.4dex for BMB 78 and BD-012971, whereas for BMB 289 the 


agreement in the metallicity is excellent (the fit to the observed 
spectrum of that star is somewhat better with the parameters of 
iRich & Origlial |200^ . However, the metallic ity of BD-012971 
was re cently also determined independently bv iRvde & SchultheisI 
( l2015h using different data to be [Fe/H] = -0.78, in good agree¬ 
ment with our result ([Fe/H] = -0.68). 

The C/O ratio found in our modelling is higher than the litera¬ 
ture value by large factors in some cases. To identify the reason for 
thi s difference, the origina l spe ctra of the comparison s tars obtained 
bv iRich & Origlial l l2005h and ICunha & Smitn ( l2006l) were kindly 
made available to us by these authors. A comparison showed that 
the spectra of BMB 78 and BD-012971 are indeed very similar to 
ours, i.e. the difference in the abundances can not be attributed to 
the different observational material. Some variatio n can be seen for 
the sp ectra of BMB 289. This star is identified bv ISoszvhski et all 
l l2013h to be a semi-regular variable with a primary pulsation pe¬ 
riod of 162.53 d and an I-band amplitude of Of 132. Most likely, 
the difference seen in fhe specfra is due to dynamic effects in the 
atmosphere of that star. The fit to the spectrum of BMB 289 is of 
lower quality t han for the other sta rs. Also, because of its late spec¬ 
tral type (M9, iBlanco et al.lll984h . this star could be considerably 
cooler than assumed in the literature. For these reasons, BMB 289 
appears to be not a good comparison star and we put less weight on 
it. 

The C/O ratios reported by iRich & Origlial (l2005h are very 
low. Using the abundance pattern given in the literature, our mod¬ 
els do not give a proper fit of the observed spectra. Note that Tts 
and log g as given in the literature were adopted. The reason for the 
differences found in the C/O ratios must lie in the different analy¬ 
sis methods. Because of the insensitivity of the OH lines available 
here to the O abundance in M-type giants (Sect. l3.3T l. we assumed a 
priori a trend of the O abundance with metallicity that was derived 
from hotter stars. This trend might not be applicable to our sample. 
However, we believe that most of the difference in the results stems 
from inconsistencies in the O abundance used for model construc¬ 
tion and for spectral synthesis, which may produce different C/O 
ratios. We should add that the measured C/O ratios in our sample 
are expected to be reliable at least in a relative sense, so that we are 
able to draw conclusions on the occurrence of internal C enrich¬ 
ment by 3DUR 

Table 1^ also collects the '^C/'^C ratios reported in the litera¬ 
ture dRich & Ori gliall?005l) and measured here. The ratio is quite 
uncertain for BMB 289, but for BD-012971 the agreement is very 
good. 

In conclusion, the agreement with literature results on the in¬ 
dividual star level is mixed. The excellent agreement that was found 
for the well-studied star Arcturus however gives confidence that the 
abundance patterns that we derived for our science targets are reli¬ 
able. As shown in Sect. 14.21 the good agreement in mean metallic¬ 
ity on the level of late-type star samples in the inner bulge regions 
gives us further confidence in our results. 


4 RESULTS AND DISCUSSION 

In most of the following subsections, the Plant stars need to be dis¬ 
cussed separately from the Spitzer stars because they stand out from 
the sample in terms of location in the bulge, but also in metallic¬ 
ity and C/O ratio. Comparisons with other bulge samples are based 
only on the Spitzer sample because of closer spatial overlap. We 
ask the reader to pay attention to which subsample is discussed in 
which place. 
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Table 4. Comparison of abundances with literature results. 


Name 

[Fe/H] 
this work 

[Fe/H] 

lit. 

C/0 

this work 

C/0 

lit. 

>2C/‘3C 

this work 

‘^C/'^C 

lit. 

Ref^ 

BMB 78 

-0.47 

-0.03 

0.46 

0.20 

9.9 


CS06 

BMB 289 

-0.05 

-0.05 

0.57 

0.28 

4.0: 


CS06 

BMB 289 

-0.21 

-0.15 

0.67 

0.12 

4.6: 

5.0 ± 1.2 

RO05 

BD-012971 

-0.68 

-0.24 

0.30 

0.15 

9.5 

7.9 ± 1.8 

RO05 

BD-012971 

-0.68 

-0.78 

0.30 


9.5 


RS15 


References: CS06 = ICunha & SmithI i2006h . RO05 = iRich & Origlial ilOOSh . RS15 = 
iRvde & Schultheisl j2015h . 


4.1 Current masses 


The current masses may tell us something about the age of the 
sample stars and the star formation history in the Galactic bulge. 
Due to stellar mass loss, the initial masses must have been higher 
than the current masses. Some of the masses listed in Table [T] are 
significantly larger than I.OMq, implying ages considerably lower 
than lOGyrs. This is in agreement with w hat was found previ¬ 
ously from studies of AGB stars in the bulge l van^oon_et alj2003l: 
[Groenewegen & Blommaertll2005l : lUttenthaler et al.ll2007h : lunii- 
nosities, pulsation periods, and dredge-up behaviour suggest that 
there are stars of masses 1.5 - 2.0Mq and intermediate age (< 
5 Gyrs) present in the bulge. A lthough not obvious in purely pho¬ 
tometric stu dies of the turn-off JZoccali et al.ll2003l : IClarkson et al.l 
l2008Ll201lh . recent investigations of m icrolensed dwarf and sub¬ 
giant stars with known metallicities by iBensbv et alj l l2013h indi¬ 
cate that a younger population is probably present in the bulge. The 
masses derived for our sample stars are in line with those results. 
An extended star formation history is required to explain the more 
massive AGB stars in the bulge. 


The four LAVs among the Spitzer sample seem to form a dis¬ 
tinct group. They have the following independent properties that 
distinguish them from the other stars: their spectra, with very strong 
TiO and VO bands, make them the coolest stars observed; their lu¬ 
minosities, if corrected for the depth scatter, would be the highest 
of the stars in the sample; they have the longest observed pulsation 
periods, from 392 to 472 days; and they have both the lowest and 
highest radial velocities of any of the stars observed, suggesting that 
their velocity dispersion is higher than that of the other stars. The 
probability to find the extreme RV values among any sub-sample 
of four stars within our Spitzer sample of 37 stars is 0.9%. The 
upper panel of Fig. [7] shows the distribution of the Spitzer sample 
stars in the RV vs. period diagram. The lower panel of that figure 
shows a band strength index (pseudo-continuum: 911 - 921 nm, 
band: 923 - 938 nm) as a funclion of period. This band slrength 
index shows a very slrong linear correlation with the determined 
temperature. The four LAVs are plotted as red symbols. The four 
LAVs not only have the highest and lowest radial velocity, they also 
have the strongest molecular bands, indicating that they form a dis¬ 
tinct group among the Spitzer sample. 


The high luminosities of these stars are consistent with their 
high derived masses in that high-mass stars evolve further up the 
AGB before they eject all their envelopes. The very low tempera¬ 
tures, especially if the stars really are massive, suggest that these 
stars are of high metallicity since the giant branch becomes cooler 
with increasing metallicity, but hotter with increasing mass. Unfor- 



0 100 200 JOO 400 500 

Period (d) 


Figure 7. Radial velocity and band strength index as a function of pulsation 
period of the Spitzer sample stars. The LAVs are shown as red symbols. The 
extreme velocities as well as the highest band strengths are found among the 
stars with the longest periods. Stars for which no period is known are set to 
P = 0. 


tunately, their near-IR spectra are too complex to derive any abun- 
danc es from them. _ 

IWood & Besselll ( Il983h observed a similar group of long pe¬ 
riod variables in the bulge and they also deduced masses of 3 - 
5Mq. They concluded that these objects had metallicities about 
2.5Zq and suggested that they could be young metal-rich stars 
ejected from star-forming regions in the Galactic plane near the 
centre of the Galaxy. Star clusters recently formed n ear the Galac¬ 
tic centre are now well known (e.g. lWang et al.l2006l) . The 3-5Mq 
stars could also be ejecte d from the central st ar cluster where young 
stars are known to exist iGenzel et alj|201oh . The high RV disper¬ 
sion obtained for the present sample of four cool, high-luminosity, 
long-period variables lends support to this suggestion. 


4.2 Metallicity 

The metallicity distribution function (MDF) of the sample stars is 
shown in the middle panel of Fig. [8] The Plant stars are included 
in this diagram. The MDF is symmetric and shows a clear peak 
between [M/H] = -0.10 and -0.20. 

The median metallicity of the Spitzer sample is [M/H] = 
-0.15, the mean is ([M/H]> = -0.14 ± 0.03, with a standard devi¬ 
ation of 0.15 ± 0.02 dex. Only few stars have super-solar metallic¬ 
ity. The Spitzer sample appears to be very homogeneous. The re¬ 
sult for the Spitzer sample may be compared to some results from 
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Figure 8. Top panel: Histogram of the C/0 ratios. Middle panel: Histogram 
of the metallicity of all sample stars with reliable metallicity determination. 
Bottom panel: C/O ratio as a function of metallicity. The Plaut stars are 
shown as red symbols, also two typical en'or bars are shown. 


the literature: iRich & Origlial ( l2005h find ([Fe/H]> = -0.19 ± 0.02 
from the analysis of high-resolution (R = 25 000) H-hand spectra 
of 14 M-type giants in Baade’s Window at {l,b) = (-1-1°,-3.9°), 
in excellent agreement with our result. The Icr dispersion in their 
sample is 0.080 ± 0.0 15 dex, somewhat narrower than in ours. 
ICunha & SmithI ( l2006h analysed high-resolution {R = 50 000) H- 
and K-hand spectra of seven bulge K- and M-type giants and derive 
([Fe/H]) = -0.15. Despite some disagreement with those works on 
indivi dual objects (Sect . lT^ . the mean metallicities of the samples 
agree. l^ch et al.l (l2007h observed 17 M giants in a field 1° south of 


Figure 9. Comparison of MDFs: this work ( s haded hist ogram); combined 
MDF from iRich & Origlial ilOOSll.l Rich et al] j2007h . and iRich et al.l j2012L 
dotted linei: and iBensbv et alJ j2013l dashed linei. 


the Galactic centre and measure ([Fe/H ]) = -0.22 ±0.03 with a Icr 
dispersion of 0.140 ±0.024 dex. Finallv. lRich et alj j2012h observed 
a total of 30 M giants in two fields 1.75° and 2.65° south of the 
Galactic centre and find in these fields ([Fe/H]) = -0.16±0.03 and 
a = 0.12±0.02, and ([Fe/H]) = -0.21 ±0.02 and tr = 0.09±0.016, 
respectively. All these results for (non-variable) K- and M-type gi¬ 
ants agree very nicely with our result for the Spitzer sample. 

A lack of metal-rich stars among M-type giants in the 
bulge was noted in previous works (e.g. lUttenthaler et al.l 12012 : 
iRich et aH 1201 2l . and references therein). This is surprising be¬ 
cause only relatively metal-rich stars would be expected to be¬ 
come cool enough to evolve to M-type giants. Figure |9] shows 
a comparison between MDFs of in total 61 bulge M-type giants 
jRich & Origliall2005l: [Rich et al.ll200l l2012h and AGB stars (27 
stars, this work), and the MD F of bulge dwarf and subgiant stars 
(58 stars. iBensbv et al.ll2013h . While the MDFs of M-giants and 
AGB stars are quite narrow and peak between [Fe/H] = -0.3 and 
-0.1, the MDF of dwarfs is much broader with less clear peaks. The 
last also contains more metal-rich stars than do the giant samples. 
A two-sided Kolmogorov-Smirnov test reveals that the probability 
that our Spitzer and the dwarf star sample are drawn from the same 
underlying population is only 0.01. For the 61 M giants and the 
dwarf sample, the mutual probability is even less than 10“®. 

It has been suggested that very high mass-loss rates at high 
metallicities cou ld remove those stars from the canonical paths o f 
stellar evolution dCohen et al.l2008l : ICastellani & Castellanilll993h . 
hence they do not reach the coolest and most advanced phases. 
The low number of super-solar metallicity stars in our sample is 
in agreement with this conjecture. 

On the other hand. iBueTl ( |20 1 3l) finds that bulge PNe descend 
mostly from super-solar metallicity stars. This result may be biased 
towards high-metallicity objects because only PNe with relatively 
high surface brightness can be studied at the distance of the bulge. 
Furthermore, metal-rich stars may be more likely to form PNe than 
metal-poor ones, e.g. due to enhanced mass loss and/or a higher 
binary fraction. 

The two Plaut stars for which a metallicity could be deter¬ 
mined have a lower metallicity than the remaining sample (one 
Spitzer sample star has the same metallicity as one of the Plaut 
stars). Though not a st rong confirmation , this fi nding is consistent 
with the conclusion of lUttenthaler et aP ( 1201 2h that AGB stars in 
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the outer bulge {b x -10°) descend mainly from the metal-poor 
bulge population, not from the metal-rich one. This conclusion was 
however solely based on the kinematical properties (radial velocity 
dispersion) of AGB stars and of the metal-poor and metal-rich pop¬ 
ulation of RGB stars in a field cent red on (I, b) = (0° , -10°) , not on 
actual metallicity determinations. lUttenthaler et all l l2012l) found 
the metal-poor population in that field to have a mean metallicity 
[M/H] = -0.57 ±0.03 with a dispersion of 0.28 ± 0.03 dex. The two 
Plant stars analysed here fall well into this distribution, which sug¬ 
gests that they actually descend from this metal-poor bulge popula¬ 
tion. Note however, that non-negligible biases could influence this 
finding: Of t he 671 long-period and semi-regular variables found by 
IPlaud jl97ll) in the Palomar -Groningen field no. 3, 2 7 were selected 
for optical spectroscopy bv lUttenthaler et al] ( l2007h . of which only 
eight have been observed with CRIRES and analysed here, four of 
them being Tc-rich. 

A negative metallicity gradient ha s previously been f ound 
along the minor axis of the Galactic bulge dZoccali et alj2008l) . Our 
complete sample is in agreement with such a gradient because the 
metal-poor Plaut stars are far from the plane. Even when only the 
Spitzer sample is considered, we find a slight negative metallicity 
gradient, but this is not very significant. 
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Figure 10. Carbon isotopic ratio vs. C/0. Filled black circles represent stars 
from this work, while filled blue triangles a re Tc-rich MS- or S-type , open 
blue triangles Tc-poor M-type giants from ISmith & Lambertl il990l) . The 
filled red diamond symbol is the Tc-rich star Plaut 3-942. The two dashed 
lines indicate the evolution of atmospheres with (C/0,'^C/'^C)=(0.30,20) 
and (0.40,10) by adding *^C to them. A typical error bar is shown for one 
star (J174917.0-293502). 



4.3 C/O and carbon isotopic ratios - Which AGB stars 
undergo 3DUP? 

All sample stars are oxygen-rich, i.e. there is no carbon star in the 
sample. Also among those sample stars whose abundances could 
not be measured, none has an enhanced strength of the CO band 
head, except for the Plaut stars with Tc. The distribution of C/O 
ratios in the sample is shown in the upper panel of Fig. [8] Our 
sample shows an even distribution of C/O ratios between 0.22 and 
0.51, again except for the Plaut stars, which have a considerably 
higher C/O ratio. This was expected because these two stars were 
found to show lines of Tc, an indicator of recent or ongoing 3DUP. 
A comparison between the spectra of the Plaut stars with those of 
similar stars from the Spitzer sample shows that the former indeed 
have much more prominent CO and CN lines. Hence, the two Plaut 
stars are a good benchmark of which C/O ratio we may expect to 
find in stars that underwent 3DUP and to check if any of the Spitzer 
sample stars underwent 3DUP. 

The mixed chemistry ob served in many Galactic bulge PNe 
jGuzman-Ramirez et alj[20llh may be understood if those objects 
descend from stars such as the Tc-rich Plaut stars. They have an 
enhanced C/O ratio (> 0.7), which would ease the formation of 
PAHs in their outflows once they become post-AGB objects or PNe. 
Our finding of a relatively high C/O ratio in two of the Plaut stars 
(C/O ^ 0.7) does not exclude the existence of an UV-irradiated 
torus around the centra l stars of Galactic bulge PNe suggested by 
iGuzman-Ramirez et al.l(l201 ih . Nevertheless, our observations sug¬ 
gest that the assumption of low C/O ratios in PNe precursors to¬ 
wards the Galactic bulge is not easily justified. _ 

The solar C/O ratio in the compilation of ICaffau et alj l l2008tl 
used here is 0.55, hence all Spitzer sample stars have an at least 
slightly sub-solar C/O ratio. The mean C/O ratio of the Spitzer stars 
is 0.38, with standard deviation of 0.10. Careful measurements of 
the C/O ratio in K-typ e bulge giants find so mewhat lower ratios 
than this. For example. iMelendez et alj 1 20081) find a m ean of 0.24 
in a sample of 19 K-type giants, and Rvde et ahl ( l2010h find a mean 
of 0.26 among 11 stars. Although the Spitzer sample has a some¬ 
what higher mean C/O ratio than these K giant samples, we do not 
think that the sample contains stars that underwent 3DUP. Rather, 


there could be systematic effects between the analysis of the AGB 
stars here and the hotter K-type giants. 

Stars that have undergone 3DUP events may be revealed by 
a plot of the carbon isotopic ratio vs. C/O, since both these ra¬ 
tios are expected to be enhanced. Such a diagram is presented in 
Fig. [To] I ncluded in that diagram are solar neighbourhood red gi¬ 
ants from ISmith & Lambert! Ill 99(il) : Filled blue triangles represent 
giants of type M or MS with Tc, open blue triangles are M-type 
giants without Tc. Tc-poor MS/S-type stars are not included here 
because they are thought to be the product of binary mass transfer. 
To guide the eye, the evolution of AGB atmospheric abundance pat¬ 
terns by adding *^C to them is indicated by the two dashed lines in 
Fig . 1101 This is a rather simplistic vi ew because it igno res possible 
effects of deep mixing on the AGB ( IBusso et alj|2010l) . The oper¬ 
ation of deep mixing can considerably mask the 3DUP evidence 
because '^C will be consumed to from '"'N. Indeed, not all of the 
Tc-rich stars in Fig.[T0|have an enhanced C/O or *^C/*'’C ratio, and 
there are Tc-poor stars that have a higher C/O ratio than some of 
the Tc-rich stars have. The disc giants discussed in the Appendix 
may give a further indication of typical abundance ratios found in 
pre-3DUP giants. However, there are a few disc stars in Fie.llOlthat 
have obviously added a considerable amount of ’^C to their atmo¬ 
spheres. The only star from our sample that is displaced from the 
rest is Plaut 3-942 (filled red diamond symbol), a star with known 
Tc content. Its '^C/'^C ratio is uncertain, but seems to be higher 
than that of most of the other sample stars. Unfortunately, it was 
not possible to determine the '^C/'^C ratio in Plaut 3-626. 

The Spitzer sample star with the highest ratios is J180308.7- 
295220 at C/O=0.51 and ‘2c/‘3C=28. Note, however, that the er¬ 
rors in these quantities are correlated, an overestimate in C/O also 
leads to an overestimate in ‘^C/'^C, see Sect. l3.5l This star does not 
show an enhanced Y abundance in its atmosphere (Table |2). Also, 
none of the cool stars for which no metallicity or C/O ratio could 
be determined has a significantly enhanced '^C/'^C. Therefore, we 
believe that none of the Spitzer sample stars has undergone a 3DUP 
event yet (although they may have undergone He-shell flashes with¬ 
out 3DUP). The same conclusion is reached from the location of the 
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Figure 11. Formation of carbon star's and 3 DUP behaviour in t he ini¬ 
tial mass vs. metallicity plane of models from iMarigo et all l2013h . Light 
coloured areas indicate a combination of mass and metallicity where the 
model star becomes carbon-rich for at least a short while on the AGB. The 
darker the colour, the lower the maximum C/O ratio reached by the model 
star. The white symbols are sample stars for which a current pulsation mass 
and a metallicity was determined: circles ai'e for Spitzer sample star's, dia¬ 
monds for Plaut stars. 


sample stars in a dust mass-loss rate {K - [22]) vs. pulsation period 
diagra m, •which allo'ws to distinguish between Miras with and with¬ 
out Tc (IUttenthaleill2013h . This also implies that C/O and 
are not useful as “evolutionary clock” for the Spitzer sample, as 
initially anticipated for the observations. 

With the sample of stars for which we have information on the 
current mass and the metallicity, we may put constraints on AGB 
evolutionary models that make predictions concerning the occur¬ 
rence of 3DUP and the formation of carbon stars. It is difficult to 
obtain the mass and metallicity of an AGB star, because it requires 
a star to pulsate in order to derive its mass, but dynamic effects 
caused by the pulsations hamper the determination of the metal¬ 
licity. Nevertheless, in Fig. m the stars are plotted in the mass - 
metallicity plane, along with p redictions from mo del tracks calcu¬ 
lated with the COLIBRI code l lMarigo et alj|2013h . The maximum 
C/O ratio reached by the model stars is colour-coded in that dia¬ 
gram: stars in the light coloured area in that plane turn into carbon 
stars at least for some time on the AGB, the darker the colour the 
lower the maximum C/O ratio. The minimum initial mass required 
for a star to experience 3DUP or to become a carbon star increases 
with increasing metallicity because the mixing is le ss efficient at 
higher metallicity. At solar metallicity, the models of iMarigo et al.l 
( I 2 OI 3 I 1 predict that a minimum initial mass of 1.90 Mq is required 
to form a carbon star, and the interpulse luminosity these stars is 
predicted to be in the range 2 OOO- 6 OOOL 0 . Our sample stars (white 
symbols) have sub-solar metallicity, but their current masses would 
be too small to undergo significant 3DUP. 

The two Plaut stars, which are at moderately low metallic¬ 
ity and clearly experienced 3DUP, are also in a region of the di¬ 
agram where the models do not predict 3DUP to happen, clearly 
in contradiction with observations. In principle, 3DUP should be 
eased in these stars because of their lower metallicity. Neverthe¬ 
less, at the metallicity of the Plaut stars, the models predict that 
a higher initial mass is required for 3DUP to occur than the Plaut 
stars have. The initial masses of the stars must have been higher 
than the current masses, but it is not expected that they have lost a 


significant amount of their mass. The largest mass loss occurs dur¬ 
ing the super-wind phase on the AGB, which our stars have not yet 
entered, except maybe the ones with the highest mass-loss rates. 
If the determined current masses are reliable and if they have not 
been lowered much from the initial value by mass loss, then the 
evolutionary models would predict 3DUP to occur at too high ini¬ 
tial masses at these moderately low metallicities. 

A number of sample stars with a high current mass but with¬ 
out measured metallicity might fall in the range where carbon star 
formation is predicted by the models. Since all of these stars are 
oxygen-rich, they are problematic for the evolutionary models, un¬ 
less they have a considerably higher metallicity than the other sam¬ 
ple stars. This may be the case for the LAVs in the Spitzer subsam¬ 
ple, see Sect. 14. II We caution that there is also the possibility that 
the linear pulsation models are not adequate to derive the current 
masses of some of the stars, or that the input stellar parameters for 
the mass determination (e.g. temperature) are in error. 

The lower panel of Fig. [ 8 ] shows the C/O ratio as a function 
of metallicity. It is clear that the Spitzer sample stars are quite ho¬ 
mogeneous in their overall composition, with metallicities mostly 
sub-solar and all C/O ratios below the solar value. The Plaut stars 
(red symbols) are clearly offset from this group. Our sample does 
however not allow to infer the evolution of the C/O ratio with metal¬ 
licity in the bulge reg ion. For mode l s and observations of the C/O 
evolution, we refer to lCescutti et alj (l2009h . 


4.4 Metal abundances 

Only the Spitzer sample stars were observed in the additional K- 
band setting; the individual metal abundances derived from them 
are discussed here. The results of these measurements are presented 
in Fig.[T2 The scatter in the abundances may be fully explained by 
measurement uncertainty (see the error bars on the left hand side 
of the panels). The abundances are compared primarily to mea¬ 
surements in similar regions of the bulge as the location of the 
Spitzer sample stars ( iBensbv et al.ll2013l ; lAlves-Brito et al.ll201(]| ; 
lOonzalez et al.ll 201 lh . Only for Al measurements f rom t he outer 
bulge (Plaut’s field, b k: - 8 °) from Ijohnson et aO ( 1201 2 h are in¬ 
cluded. The Spitzer sample is too small and spans too small a range 
in metallicity to probe for any trends of the elements with metallic¬ 
ity. However, the abundances are suited to check if the AGB stars 
in the present sample descend from the same populations that can 
be found among less evolved bulge stars. A look at Fig.lJ^can an¬ 
swer this question with “yes”: The relative abundances of Al, Si, 
Ti, and Y agree well with the range of abundances that are seen in 
other bulge stars at these metallicities. The fact that Y, one of the 
light s-process elements, is not significantly enhanced also argues 
against 3DUP being active in the Spitzer sample stars for which 
the Y abundance could be derived. Also, no correlation between 
[Y/Fe] and C/O is found in our data. C and Y are thus probably not 
produced in situ and the abundances likely reflect the primordial 
values. 


4.5 Radial velocities 

Here only the heliocentric radial velocities of the 37 Spitzer sample 
stars (column 2 in Table are discussed. For t he radial velocity 
disper sion of AGB stars in the outer bulge, see lUttenthaler et ^ 
1 I 2 OII Section 4.7.3). 

The radial velocities of the Spitzer stars have a mean of 
(RV) = 16.1 ± 15.8 km s“' and a dispersion of ctrv = 95.8 ± 
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Figure 12. Abundances of [Al/Fe], [Si/Fe], [Ti/Fe], and [Y/Fe] as a func¬ 
tion of [Fe/H] (from top to bottom) in samples of Galactic bulge stars. 
Filled black circles are data points from this work , black plus symbols are 
dwarf and subgiant stai's fromiBen sbv et al ]i2Q13h . blue open diamonds ai'e 
from Baade’s field fro m lAlves-Brit^^ aHiMS, green open circles from 
iGonzalez et 7i H20iiL only Baade's field), and red crosses are from Plant’s 
field in the outer bulge from ljohnson et al] i2Q12h . The big symbols on the 
left hand side of each panel show typical error bars. 


11.3 km s ‘. However, it appears that the RV dispersion depends on 

the latitude of the fields (Fig. [TSl l: while in the field at fc-3.8° the 

dispersion of the 17 stars is only 65.9 ± 11.6 km s“*, the fields at ~ 
±1° (15 stars) have a combined dispersion of 115.0 ± 21.7kms“'. 
It has been shown that the metal-rich stellar population of the 
Galactic bulge has a radial velocity dispersion that significantly in¬ 
creases when approaching the plane, whereas the metal-poor pop- 
ulation has a con s tant q -r v of 90 - 100 km s ~‘ thr o ughout the bulge 
I Babusiaux et al. I 20 id : lUttenthaler et al.l l2012l : iBabusiaux et al.l 
l2014h . This increasing radial velocity dispersion is a result of 
the action of the Galactic bar. We thus suggest that the Spitzer 
sample also follows a bar-like kinematic pattern, and that it be¬ 
longs to the metal-rich population. The metallicities measured in 
the present paper are in agreement with this picture. Note that 
Mira stars in the bulge are actua lly found to trace a tilted bar 
dGroenewegen & BlommaerJ20051) . As shown above that the Plant 
stars are more metal-poor. The radial velocity dispersion of AGB 
stars in the outer bulge suggests that they descend from the metal- 
poor, spheroidal component of the bulge. 


5 CONCLUSIONS 

We present an analysis of high-resolution near-IR spectra of a sam¬ 
ple of Galactic bulge AGB stars. The sample consists of a larger 



Figure 13. Radial velocity vs. absolute value of Galactic latitude for the 
Spitzer sample stars. Red symbols represent stars located north of the 
Galactic plane, at positive latitude. The large symbols next to the three lati¬ 
tude groups show the mean and standard deviation of the radial velocity of 
the respective group. 


sub-sample from the inner and intermediate bulge (Spitzer sample) 
and a smaller sample from the outer bulge (Plant stars). Using state- 
of-the-art hydrostatic model atmospheres for cool stars and the best 
available atomic and molecular line lists for the observed wave¬ 
length range, we are able to successfully model the spectra of stars 
down to ^eff ~ 3100 K. General metallicities [Fe/H], the carbon-to- 
oxygen ratio C/0, carbon isotopic ratio '^C/'-’C, the abundance of 
Si, Al, Ti, and Y, and heliocentric radial velocities are determined. 
The abundance analysis relies on the assumption of an [a/Fe] trend 
with [Fe/H] because of the low sensitivity of OH lines to the O- 
abundance in cool, M-type giants. We caution that calculation of 
both the model atmosphere and the spectral synthesis with consis¬ 
tent abundances is required in such cool stars for a reliable mod¬ 
elling. Current stellar masses are also derived from linear pulsation 
models. 

This paper touches on the problem of the apparent lack of 
intrinsic carbon stars in the Galactic bulge. All sample stars are 
oxygen-rich, there is no carbon star among them. The C/0 and car¬ 
bon isotopic ratios measured in the sample stars suggest that the 
Spitzer sample stars did not undergo the third dredge-up. They ap¬ 
pear to have a composition typical for post-first dredge-up red gi¬ 
ants. Furthermore, no increased CO band head strength was found 
among the stars for which we could not measure the overall compo¬ 
sition, excluding the known technetium-rich stars. Note, however, 
that the effect of 3DUP on the carbon abundance may be consid¬ 
erably masked by the operation of deep mixing on the AGB. The 
conclusion of a lack of 3DUP in these stars is supported by the fact 
that Y is not enhanced in the stars with observations of a line of that 
s-element, and the posit ion of the Spitzer stars in a dust mass-loss 
rate vs. period diagram ( IUttenthaleill2013h . On the other hand, two 
of the Plant stars have a clearly enhanced C/0 ratio, as expected 
from the fact that Tc lines have been found earlier in their optical 
spectra. 

Howeve r, det ailed AGB evolutionary models, e.g. 
iMarigo et alj i2013h . are in conflict with our measurements 
of metallicities, masses, and 3DUP occurrence in the sample stars. 
Those models predict no 3DUP for stars with initial masses as low 
as the current masses of some of the Tc-rich Plant stars, but do 
predict 3DUP and carbon star formation for the more massive stars 
in our sample. Though it is impossible to determine the metallicity 
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and C/O ratio for the coolest and putatively most massive stars in 
the sample, it is clear that they are not carbon-rich. Assuming a 
similar metallicity distribution as the other Spitzer stars, one would 
expect from the models that they evolve into carbon stars. The 
reason for this could be that the masses estimated with our linear 
pulsation models are not accurate enough because of uncertainties 
on the input parameters and/or invalid assumptions for the models. 
We suggest that 3DUP is found in the Plaut stars but not in the 
Spitzer stars because of their lower metallicity, see below. On the 
other hand, the C-star lifetime is very short at high metallicity, 
and also a high He abun dance could prevent 3DUP and C-star 
formation ( lKarakasll2014h . 

It has been suggested that a high oxygen abundance could 
preve nt the formation of carbon stars in the Galactic bulge Jpeastl 
l2007h . However, the current consensus is that the O-abundance 
trend in bulge stars is very similar to that of local thick disc stars 
jMelendez et aljEoOSi : iRvde et aljEoiOl : ISensbv et alj|2013h . We 
suggest that at the slightly subsolar metallicity of the Spitzer stars, 
carbon star formation is suppressed because 3DUP is suppressed al¬ 
together, either because of the high metallicity or possibly because 
of a high He abundance. Hence, an enhanced oxygen abundance 
is not required to explain the lack of carbon stars in the Galactic 
bulge. 3DUP does occur at the lower metallicity of the Plaut stars 
in the outer bulge, but there may not be sufficient 3DUP events to 
turn them into C-rich stars, and/or because of the action of deep 
mixing. From the masses that we derive for these stars one might 
expect C-star formation, but note the uncertainties in the mass de¬ 
terminations. 


The Spitzer sample stars have on average slightly sub-solar 
metallicity but there are few that have super-solar metallicity. This 
is in agreement with the results of previous observations of cool 
giant samples in the Galactic bulge, though the comparison on the 
level of individual stars is more uncertain. The lack of super-solar 
metalli city stars could be caused by metallicity-dependent stellar 
winds dCohen et Ii]l2008b . 


The Plaut stars are found to be more meta l -poor, in agree¬ 
ment with earlier suggestions bv luttenthaler et al] 1 I 2 OI 2 I 1 that were 
based on kinematical arguments. Thus, we suggest that the AGB 
stars in the outer Galactic bulge mainly descend from the metal- 
poor bulge population ([Fe/H] ~ -0.6), not from the metal-rich 
one ([Fe/H] — hO.3), although about 30% of the RGB stars in a 
field at (I, b) = (0°, -10 °) belong to the metal-rich bulge population 
JUttentha ler et al.l 2012h . The metallicity inferred for Galactic bulge 
PNe I Buelfco" ) is higher than what we find among AGB stars. It 
could be that metal-rich stars are more likely to form PNe, which 
would be in agreement with the call for metallicity-dependent stel¬ 
lar winds. 


The abundances of Al, Si, Ti, and Y measured for some of 
the Spitzer sample stars suggest that they descend from the bulk 
of the population in their bulge region. The radial velocity disper¬ 
sion of the Spitzer stars suggests that they follow a bar-like kine¬ 
matic pattern, characteristic for the metal-rich bulge population. 
The Plaut stars, on the other hand, seem to belong to the metal- 
poor, spheroidal bulge component. 

Finally, the masses of the sample stars suggest that the Galac¬ 
tic bulge is not exclusively old, a younger and more massive popu¬ 
lation (~ 1.3Mq) seems to be present. This would be in agreement 
with age estimates of metal-rich dwarf stars in the Galactic bulge 
jBensbv et al.ll2013h . 
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APPENDIX A: ABUNDANCES IN FIVE FIELD RED 
GIANTS 

Observations of five red giant stars in the field have been obtained 
in the same CRIRES run. These spectra were also analysed and the 
results are reported here for the sake of completeness. 

The determination of the effective surface temperature and 
gravity of the field stars relied purely on literature val ues and near- 
IR pho tometry. Temperatures of the stars are given bv lWright et al.l 
(l2003h . which are based on spectral type determinations. As a 


second estimate, we 
iLebzelter et alj (12012 

applied the T^tf - {J - K)o calibration of 
) to the (J - K'fa colour that was derived from 

2MASS photometry ( 

Skrutskie et alj2006h and the interstellar red- 

dening mao ofischlegel et al. ( 1998(1. A compromise value between 


these two estimates was adopted for the model atmospheres. The 
surface gravit y log q was estimated also with the log g - {J - K)o 
calibration of iLebzelter et alj (l2012h . The stellar parameters of the 
field giants are summarised in Table lAll 

The CRIRES spectra of the field stars, which are of high qual¬ 
ity, were analysed in the same way as the bulge AGB stars in the 
main part of the paper, except that for the ff-abundance a more disk¬ 
like relation was adopted, with [tr/Re] = -1-0.4 up to [Re/H] = -1.0, 
followed by a decline of 0.4dex/dex for higher metallicities. The 
results of the analysis are summarised in Table lA2l All five stars 
have close t o solar metallicity, and slightly sub-solar C/O ratios (on 
the scale of lCafFau~et '^l200i . 
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Table Al. Stellar parameters of the field red giants. Meaning of th e columns: object na me; range of spectral types found in VizieR; J, H, and K 
magnitudes from 2MASS; interstellar e xtinction in the V-band from ISchlegel et al]il998h : de-redden ed J-K colour; effe ctive temperature derived 
from the T^ff - (J — K)o calibration of iLebzelter et “51 i2Q12l L12); effective temperature listed in lWright et alj i2QQ:l W0 3): adopted effective 
temperature; adopted log g of the model atmosphere, estimated with the log g - {J - K)o calibration of iLebzelter et a 


Object name 

Sp. type 

J 

(mag) 

H 

(mag) 

K 

(mag) 

Av 

(mag) 


Teff 

(K, LI2) 

^cfF 

(K, W03) 

T'eff 

(K, adopted) 

log 9 
(cgs) 

HD 136546 

K5III 

4.445 

3.515 

3.313 

0.730 

1.021 

3785 

3950 

3900 

1.20 

HD 140672 

K5-M11I1 

4.207 

3.251 

3.021 

0.947 

1.042 

3740 

3720 

3750 

1.10 

HD141189 

MO-2 111 

4.542 

3.551 

3.252 

0.407 

1.228 

3295 

3620 

3400 

0.38 

HD141311 

K5-M1 III 

4.849 

3.862 

3.701 

0.776 

1.030 

3766 

3720 

3750 

1.10 

HD141938 

MO-2 III 

4.314 

3.367 

3.133 

0.574 

1.094 

3624 

3720 

3650 

0.90 


Table A2. Radial velocities and abundances measured for the field red giants. Meaning of the columns: ob¬ 
ject name; heliocentric radial velocity; metallicity; C/0 ratio; carbon isotopic ratio, and relative abundances 
ofAl, Si, Ti, andY. 


Object name 

helio 
(km s“‘) 

[Fe/H] 

C/0 


[Al/Fe] 

[Si/Fe] 

[Ti/Fe] 

[Y/Fe] 

HD 136546 

-9.7 

-)-0.05 

0.38 

15.4 

+0.11 

+0.05 

+0.04 

-0.18 

HD 140672 

-27.3 

+0.04 

0.43 

15.6 

+0.20 

-0.02 

+0.12 

-0.20 

HD141189 

-97.3 

-0.09 

0.28 

14.3 

+0.03 

-0.07 

+0.14 

-0.28 

HD141311 

64.6 

-0.23 

0.50 

12.3 

+0.15 

+0.06 

+0.27 

-0.17 

HD141938 

-38.7 

-0.05 

0.30 

14.4 

+0.05 

+0.08 

+0.02 

-0.08 
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